
  #include "stdio.h"
  #include "4trmul.c"
  main()
  { int i,j;
    static double a[4][5]={ {1.0,3.0,-2.0,0.0,4.0},
                            {-2.0,-1.0,5.0,-7.0,2.0},
                            {0.0,8.0,4.0,1.0,-5.0},
                            {3.0,-3.0,2.0,-4.0,1.0}};
    static double c[4][3],b[5][3]={ {4.0,5.0,-1.0},
                      {2.0,-2.0,6.0},{7.0,8.0,1.0},
                      {0.0,3.0,-5.0},{9.0,8.0,-6.0}};
    trmul(a,b,4,5,3,c);
    for (i=0; i<=3; i++)
      { for (j=0; j<=2; j++)
          printf("%13.6e  ",c[i][j]);
        printf("\n");
      }
    printf("\n");
  }

